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I.  INTRODUCTION 


A  serious  problem  In  the  theory  of  the  turbulent  flow  of  two-phase 
systems  is  that  while  very  often  one  of  the  phases  is  a  bona  fide  fluid, 
the  other  one,  if  isolated,  would  have  a  collision  mean  free  path  so 
large  that  the  usual  continuum  approximation  cannot  be  applied  to  its 
description.  Examples  of  such  a  situation  include  flows  in  pulverized 
coal  gasifiers,  certain  problems  occurring  in  meteorology,  such  as  flows 
in  clouds  and  many  others.  In  this  paper  we  suggest  a  description  of 
such  flows  by  means  of  a  mixed  system  of  equations:  one  describes  the 
flow  of  the  possibly  dilute  (generally  particulate)  phase  by  means  of  a 
transport  equation  of  the  Boltzmann  type,  while  using  the  Navier-Stokes 
(NS)  equations  for  the  fluid  phase.  Generally,  however,  one  cannot  use 
either  the  Boltzmann  equation  or  NS  in  their  "textbook  form".  First  of 
all,  the  interaction  between  the  phases  is  a  dissipative  one;  hence  the 
volume  element  of  the  phase  space  is  not  constant  in  time  and  the 
Boltzmann  equation  has  to  be  appropriately  modified. 

Second,  we  use  a  statistical  description  of  the  turbulent  flow  of 

the  fluid  phase1  and  thus,  smooth  solutions  of  NS  are  used  only  to 
describe  the  average  flow  whereas  the  probability  distribution  of 
velocity  fluctuations  is  given  by  an  expression  derived  in  a  previous 

,  2 
paper  of  ours. 

The  paper  is  organized  as  follows.  In  the  next  section  we  derive 
the  modifications  necessary  to  transport  equations  of  the  Boltzmann  type 
in  the  presence  of  dissipative  forces.  We  derive  an  explicit  expression 
for  the  extra  term  appearing  in  the  transport  equation  in  terms  of  the 
force  acting  upon  a  single  particle.  In  Sec.  3  we  briefly  review  the 
formalism  developed  in  Ref.  2  and  apply  it  to  the  coupled  two-phase 
system.  In  particular,  we  derive  moment  equations  for  the  particle 
phase  in  the  leading  approximation  of  a  Chapman-Enskog  expansion.  The 
results  are  discussed  in  Sec.  4;  in  particular,  we  point,  out  possible 
directions  in  which  our  results  can  be  generalized. 

II.  PARTICLE  TRANSPORT  IN  A  DISSIPATIVE  MEDIUM 

Consider  a  system  of  particles  with  an  ensemble  average  of  the 

phase  space  density  given  by  fix1,  v1 ) ,  where  x1  and  v1  stand  for  the 
components  of  the  particle  position  and  velocity  vector,  respectively. 
The  integral  of  f  over  the  entire  phase  space  gives  the  total  number  of 
particles.  The  particles  are  acted  upon  by  "slowly  varying"  forces 
exerted  by  the  carrying  medium:  the  interaction  between  the  particles 


V  DeUominicis  and  L.  Peliti,  "Field  Theory,  Renormalization  and 
Critical  Dynamics  Above  :  Helium,  Antiferromagnets  and  Liquid-Gas 

Systems,"  Phvs .  Rev. .  Vol .  B18,  pp .  353-376,  1978,  and  references 
quoted  there. 

2 

G.  Domokos ,  S.  Roves i-I)omokos  and  C.  K.  Zoltani,  "Random  Systems, 
Turbulence  and  Disordering  Fields,"  Phvs.  Rev.  A.  submitted  for 
publ i cation. 


is  described  by  means  of  a  collision  functional.  In  the  Boltzmann 

3 

approximation  to  the  transport  problem,  cf.,  e.g.,  Huang  ,  Chapter  5, 
the  change  in  time  of  the  phase  space  density  Is  given  by  the  equation: 

^-(fd*xd*v)  =  F[f ]d*xd*v  ,  (2.1) 

where  F  stands  for  the  collision  functional  and  the  time  derivative  is 
the  total  derivative.  In  this  paper  we  assume  for  the  sake  of  sim¬ 
plicity  that  the  phases  are  non-reactive.  Hence,  in  particular,  the 
particle  mass  (or,  in  an  inhomogeneous  particle  phase,  the  average 
particle  mass),  m,  stays  constant.  In  the  present  work,  we  need  not 
specify  a  detailed  form  of  the  collision  functional,  F [ f ] .  In  general, 

its  precise  form  appears  to  be  of  little  immediate  relevance,  as  long  a* 
the  average  momentum  transfer  in  a  binary  collision  is  nonzero.  The 
"textbook  form”  of  the  Boltzmann  equation  is  arrived  at  by  using  the 
constancy  in  time  of  the  volume  element  in  phase  space  (Liouville's 
theorem)  and  dividing  out  by  the  phase  space  volume  element.  In  the 
problem  we  are  interested  in,  Liouville's  theorem  is  not  valid:  the 
force  by  means  of  which  a  carrying  fluid  acts  upon  the  particles  is  a 
velocity  dependent  one.  Let  us  assume  that  the  equations  of  motion  are 
of  the  form: 


dx  i 

cf  r  =  V  ' 


(22) 

With  this,  one  readily  finds: 

d  a  f  1 

J-(d*xd*v)  =  pL.)/m  .  (2.3) 

av1 

Hence,  the  correct  form  of  the  transport  equation  describing  the  time 
change  of  the  phase  space  density  becomes: 


af  i  df 
_  +  v  — + 
at.  i 

ax 


*  t <*4>  -  pm 

m.i  m  „  l 

a  v  av 


(2.4) 


If  all  forces  were  conservative,  the  last,  term  on  the  l.h.s.  of  fchj. 
(2.4)  would  vanish. 

We  model  the  interphase  interaction  as  follows.  Assuming  that,  tin- 
particle  phase  consists  of  particles  which  have  random  shapes  and  are 
randomly  oriented,  it  is  clear  from  Galilean  invariance  that  the  fore: 
which  acts  upon  a  particle  of  a  given  velocity  must  be  proportional  to 
the  relative  velocity  of  the  particle  and  the  surrounding  fluid  eleir<">' 


K.  Huang, 


lanlcs .  Wiley,  New  York,  1978. 
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We  write  the  force  in  the  form: 


f1  =  671  u  a  pf( u1  -  v*)D(R  )  . 


where  the  components  of  the  fluid  velocity  are  denoted  by  u  , 


(2.5) 
Pf  stands 


for  the  mass  density  of  the  fluid,  P  is  an  effective  kinematic  viscosity 
acting  between  the  fluid  and  the  particle.  The  dimensionless  function 
D(R  )  describes  the  deviation  of  the  force  from  Stokes’  law.  It  depends 


on  the  dimensionless  "particle  Reynolds  number”, 


where 


a  is  the  average  radius  of  a  particle.  The  normalization  has  been 
chosen  in  such  a  way  that  D(0)  =1.  At  this  stage  we  do  not  specify  the 
function  D  any  further.  For  high  values  of  R^ ,  there  exist  no  good 

theoretical  expressions  for  it;  probably,  one  has  to  be  satisfied  with  a 
semiempirical  relation.  Corresponding  to  a  force  given  by  Eq.  (2.5), 
there  is  a  force  acting  upon  the  fluid,  its  density  being  given  by: 


y1  =  -  /dJv  f1  .  (2.6) 

Using  Eq.  (2.5),  a  straightforward  computation  leads  to  the  expression 
of  the  divergence  appearing  in  the  transport  Eq.  (2.4).  We  find: 


ef 


i 


=  -6n  u  a  pf(3D  +  RpD ' ) 


(2.7) 


where  the  prime  denotes  the  derivative  of  D  with  respect  its  argument. 

We  see  that  even  if  the  drag  is  assumed  to  be  obeying  Stokes'  law,  i.e., 
R  is  small,  there  is  a  finite  correction  to  the  transport  equation 

arising  from  the  shrinkage  of  the  phase  space. 

A  further  remark  is  in  order  at  this  point.  As  discussed  in  the 
Introduction,  we  are  proposing  to  apply  a  statistical  theory  of  the 

turbulent  flow  of  the  fluid  phase.  As  a  consequence,  the  components,  u1, 
of  the  fluid  velocity  are  to  be  treated  as  random  variables.  Thus,  the 
transport  equation  has  to  be  averaged  over  the  probability  distribution 
P  i 

of  u  . 


III.  STATISTICAL  THEORY.  MOMENT  EQUATIONS 

We  start  by  reviewing  the  statistical  theory  of  a  chaotic  process 
(turbulence  in  particular)  as  formulated  in  Ref.  2.  We  need  a  formula¬ 
tion  more  general  than  the  one  just  suitable  for  the  investigation  of  a  1 

single  phase  system.  I 

Consider  a  dynamical  system  which  exhibits  chaotic  behavior.  We  ; 

denote  the  collection  of  dynamical  variables  describing  the  system  by  X. 

In  practically  important  cases  X  is  an  element  of  a  vector  bundle  over 

3  1 

space-time,  i.e.,  R  x  R  .  We  assume  that  the  dynamical  equations 

satisfied  by  X  are  of  the  form,  ‘ 

I 


where  if  is  a  map  of  the  fibration  onto  itself.  Following  the  ideas  of 
Martin  et  al . .  and  DeDominicis  et  al . ,  Ref.  1,  it  was  found  in  Ref.  2 
that  the  probability  measure  over  the  (generally  infinite  dimensional) 
phase  space  characterizing  chaotic  behavior  is  governed  by  a  probability 
measure  proportional  to: 

DX  Det  (<3  .  -  «=)  exp  -  <3  .  X  -  If  \d_  -  ?>  ,  (3.2) 

t  oX  t  '  t 

where  DX  stands  for  a  volume  element  over  the  bundle  and  <  |  >  is  a 
scalar  product  over  it.  The  determinant  has  to  be  given  an  appropriate 

4 

meaning,  typically  in  terms  of  an  integral  over  Grassmann  variables. 

Our  coupled  system  of  equations,  consisting  of  Eqs.  (2.4)  and  (2.5) 
and  NS  with  an  external  force  density  given  by  Eq.  (2.6),  does  not  fit 
directly  into  this  formalism.  However,  a  Chapman-Enskog  expansion  of 
the  transport  equation  ( c f . ,  Huang,  loc .  ci t .  )  leads  to  a  system  of 
equations  which  is  of  the  form  Eq.  (3.1);  moreover,  it  also  provides  a 
practical  approach  to  the  solution.  In  this  work,  we  restrict  ourselves 
to  the  leading  approximation  of  the  Chapman-Enskog  expansion.  To  that 
order  (local  equilibrium),  the  physics  is  entirely  governed  by  the 
conservation  laws  of  particle  number,  energy  and  momentum.  Due  to  the 
fact  that  the  transport  equation  is  not  entirely  of  the  standard  form, 
we  also  have  to  rederlve  the  conservation  laws. 

Let  us  introduce  some  notation.  The  number  density  of  particles  is 
denoted  by  n,  and  if  A  is  any  function  on  the  phase  space,  its  velocity 
average  is  denoted  by  angular  brackets: 

n  =  /d*v  f,  <  A  >  =  [/d*v  f  A]/n  .  (3.3) 

Let  Y  denote  any  quantity  conserved  in  the  interpartiele  collision, 

1 .  e .  , 


Jd3  v  Y  F[f]  =  0,  (3.4) 

of.,  Huang,  loc.  cit. .  Then,  in  view  of  eqs.  (2.4)  and  (3.4)  we  have 


/d*  V  Y  +  V1  +  (-£- )  — r  +■  (^-r)/m)  f  =  0. 

at  l  mi  l 

3x  3v  3v 


(3.5) 


On  performing  some  integrations  by  parts  in  order  to  bring  the  deriva¬ 
tives  over  to  Y,  we  obtain  the  conservation  law  in  the  form- 


-jr-tn  <Y> )  +  — .(n  <v  Y> )  -  n  <v  — r>  -  (— )<f  r>  =  0.  (3.6) 

at  .  i  .  i  m  ^  l 

ax  ax  av 


P.  Ramond,  Field  Theory.  A  Modern  Primer.  Benjamin,  New  York,  1981 
Chapter  VII, 


'immgrwnBnjf mnamwrtnmrvn^mm  W  W  BE  W  W  >.■  H»  itf  vj»  g*  gr  m  v  ^  yj  > 


The  remarkable  fact  to  be  noticed  is  that,  although,  In  the  course 
of  integrating  by  parts  the  third  term  of  Eq.  (3.5),  one  would  pick  up  a 

term  proportional  to  the  derivative  of  f*  with  respect  to  the  velocity, 
that  term  cancels  exactly  with  the  last  term.  Thus,  the  expression  of 
the  consevation  laws  is  unchanged  whether  or  not  the  force  is  velocity 
dependent.  (The  standard  textbook  derivations  of  the  Boltzmann  equation 
and  the  equations  of  hydrodynamics  therefrom,  typically  assume  that  the 
forces  are  conservative.)  It  is  almost  needless  to  remark,  however, 
that  beyond  the  leading  order  of  the  Chapman-Enskog  expansion,  the 
expression  of  the  pressure  tensor  will  be  modified  due  to  the  velocity 
dependence  of  the  forces. 

We  are  now  ready  to  cast  our  equations  into  a  form  given  by  Eq. 
(3.1).  In  the  approximation  we  are  dealing  with,  and  on  denoting, 


■iwmir.f.' 


1 

>' V  >  =  w 


i 


the  dynamical  equations  are: 


Dnw  .1.  -i 

rn  -rrr —  +3Q-n<f>=0, 


Dt 


(3.7a) 


rx  v.J  .  -1  »  _i  .  ,n  »  ,*i. 


(3  +  u  d  )u  +  p  a  p  -  u  v * u  +  ( — )  <f  >  =  0.  (3.7b) 

t  r  I  P j. 


(The  somewhat  asymmetric  appearance  of  the  coupling  term  between  the  two 
equations  in  Eq.  (3.7)  is  due  to  the  fact  that  originally  we  defined  the 
interphase  coupling  in  terms  of  a  force  acting  upon  a  single  particle.) 

In  these  equations,  a^  =  d/dy}  an  Q  stands  for  the  pressure  exerted  by 

the  particuiate  phase  upon  itself.  In  most  cases,  unless  the  particle 
density  is  very  high,  one  may  set  Q  =  0.  and  we  shall  do  so  in  what 
follows.  The  symbol  D/Dt  stands  for  the  Galilei  invariant  time 

derivative,  D/Dt  =  a^  +  w1^;  on  writing  down  the  NS  equation  for  the 

fluid  phase,  we  implicitly  assumed  incompressibility.  This  cannot  be 
done  for  the  particle  phase:  the  equation  of  continuity  has  to  be  taken 
into  account  as  a  constraint. 

In  order  to  take  the  equation  of  continuity  into  account,  we  notice 
that  the  equation  of  continuity  for  the  particle  phase  reads  as  follows: 


1 


V 

wl 


afcn  +  a.fnv1)  =  0  .  (3.8) 

In  the  statistical  theory,  we  will  have  to  integrate  over  n  and  v1  with 

a  volume  element  proportional  to  Dn  Dv* .  In  the  standard  fashion,  (cf. 
Ramond,  loc .  cit . )  this  volume  element  will  be  replaced  by 

OoDr  Dn  Dv1  6(8^  +  (8^vr)n  +  v^^n]  exp  Jd’xdt  ,  (3.9) 


1 


where  a  and  r  are  Fadeev  Popov  glicSt  fields,  hence, 
algebra  valued.  The  ghost  Lagranglan,  L  ,  is  given 


they  ui e  Grassmann 
by  the  functional 
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derivative  of  the  continuity  equation  with  respect  to  n,  sandwiched 
between  the  ghost  fields: 


L  =  0(<3  +  (3  vr)  +  vrd  )T  . 

g  t  r  r 


(3.10) 


The  presence  of  a  delta- functional ,  Of...]  ,  in  the  volume  element  Eq. 
(3.9),  makes  computations  rather  awkward.  The  standard  trick  familiar 
from  gauge  theories,  namely  a  Gaussian  averaging  over  a  family  of  gauge 
conditions,  is  not  applicable  here,  because  Eq.  (3.8)  is  not  a  gauge  con¬ 
dition.  However,  one  may  replace  the  delta-functional  by  a  Gaussian,  viz . 

+  <3 ^ ( n vr )  ]  exp  -A  2  Jd*xdt  [3^n  +  d^fnv1")]2  ,  (3.11) 

where  A  is  a  parameter,  to  be  let  go  to  zero  at  the  end  of  the  calcu¬ 
lation.  Furthermore,  the  Gaussian  functional  in  Eq.  (3.11)  needs  not  be 
normalized,  since  the  physically  interesting  quantities  are  expressed  in 
the  form  of  cumulants:  their  calculation  is  not  affected  by  field- 
independent  factors  in  the  probability  measure.  Combining  now  these 
results,  we  obtain  the  constrained  functional  volume  element,  D* ,  for 
the  stochastic  (infinite  dimensional  Wiener)  integration  over  the 
dynamical  variables  of  the  particle  phase: 


I)4>  =  DoDrnv  fin  exp  T  , 


(3.12) 


where 


r  Jd’xdt{  -  A  2[0^n  +  8r(nvr)  |2  +  *  vrd^  +  e^v^r)  .  (3.13) 

Given  the  fact,  that  the  fluid  phase  has  been  treated  as  an  incom¬ 
pressible  one,  the  equation  of  continuity  needs  not  be  taken  into 
account  by  means  of  the  Fadeev -Popov  procedure.  All  one  needs  is  to 
project  out.  the  transverse  component  of  t he?  hydrodynami cal  velocity 

field,  w* ,  which  can  be  done  in  a  well-known  way. 

In  order  to  complete  the  system  Eq .  (3.7).  one  has  to  compute  the 

quantity  <f  >.  Given  a  specific  model  for  the  drag  function  in  F.q . 

(2.5),  this  may  or  may  not  be  computed  in  a  closed  form;  nevertheless, 
the  requirement  of  Galilean  invariance  tells  us  that  the  final  expres 
si  on  can  again  be  cast  in  a  form  analogous  to  Eq.  (2.5),  with  the 

hydrodynamic  velocity,  w: ,  replacing  the  particle  velocity.  (Of  course, 
the  explicit  expression  of  a  dimensionless  drag  function  in  terms  of 

| u 1  w1  j  will  be  different  from  the  originally  assumed  form  occurring 

i n  Eq .  (2.5).) 

Following  now  the  general  procedure  outlined  in  Ref.  2,  we  have  to 
construct  a  probability  measure  in  the  six-dimensional  vector  space 

formed  by  the  components  of  u*  and  w'  Given  the  fact  that  the  system 
Kq.  (3.7)  consists  of  a  set  of  coupled  equal  ions,  a  consistent  procedure 
is  to  explicitly  construct  the  functional  IF  in  Eq.  (3.1).  To  this  end. 
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we  define  the  components  of  the  vector  X  (1  <  A  <  R)  as  follows. 

First,  we  define: 

XA  =  uA  (1  <  A  <  3).  XA  =  wA”3  (3  <  A  <  6)  .  (3.14) 

In  a  similar  fashion,  we  define  embeddings  of  the  pressure  gradient  and 

the  viscous  damping  terms  into  the  six  dimensional  velocity  space.  This 

is  straightforward:  we  define  a  diagonal  tensor  projecting  onto  the 

fluid  phase:  PAI  =■  6AI>  if  1  <  A,  B  <  3,  P>r^  -  0  otherwise.  With  this, 
AB  AB  *  *  AB 

we  define: 


V  -  o  v  v2  P 

ab  pf  ; 


PA  -  0AP  (1<  A  <  3),  PA  -  0  otherwise. 


(3.15) 


In  order  to  proceed  further,  we  have  to  define  a  constant  tensor,  M, 
with  components: 


MAB  =  Vr  3  (A  <  ">•  MAB  -  V3-B  (A  >  B)' 


M  ,  0  otherwise. 

All 


(3.16) 


In  t  iirins  of  these  quantities,  we  can  define  the  six-vector  of  relative 


velocity,  Y  ,  as  follows: 


A  A  A  R 
Y  --  X  MpX 


(3.17) 


One  readily  verifies,  a  -  2 ( a J  -  v 1 ) ( u .  v.).  Thus,  the  particle 

Reynolds  number  can  be  written  in  terms  of  six  vectors  as 
1  /2 

K  (a/2  o)|y|.  The  next  step  is  to  express  the  interphase  coupling 

term;  this  is  done  by  constructing  the  vector: 


<fA>  (1  <  A  <  3).  d'1  :  <fA  3>  (4  <  A  <  6 ) 


(3.18) 


Finally,  the  embedding  of  the  gradient  operator  into  the  six  dimensional 
space  is  given  by : 


da  aA  (l  <  a  <  3),  n  -  aA  3  (4  <  a  <  6) 


(3.19) 


Willi  these  quantities,  the  system  Hq  (3.7)  ean  he  written  as  follows: 

<VH '  *\  *  '  'r*\  - 


(Pf  nm )  P(,  )X 


(V:  .  PA  .  v^xc  .  dA  =  0 


(3.20) 


This  form  allows  us  to  road  off  the  form  of  the  functional  V  and 
thus  constuct  the  probability  measure  F.q .  (3.2).  The  volume  element  is 

the  product,  1)$  Dw1  . 


IV.  DISCUSSION 

Our  main  purpose  in  this  work  has  been  to  establish  a  formalism  for 
the  dynamical  description  of  the  turbulent  flow  of  a  two- phase  system  in 
which  the  particle  phase  is  of  an  arbitrary  concentration.  The  applica¬ 
tion  of  the  Chapman-Enskog  expansion  to  the-  Boltzmann  equation,  on  the 
one  hand,  provides  a  powerful  practical  approach  to  the  solution,  on  the 
other  hand,  at  any  order  of  the  expansion  it  casts  the  equations  into  a 
form  which  can  be  handled  by  means  of  standard  techniques  of  quantum 
field  theory,  including  resummation  techniques  of  perturbation  expan¬ 
sions  such  as  the  use  of  the  renormalization  group,  etc.  In  order  to 
apply  the  present  formalism  to  realistic  cases,  such  as  flows  in  a  pipe, 
one  has  to  find  an  average  flow  around  which  one  carries  out  a  Reynolds 
decomposition  of  the  fields.  If  one  restricts  oneself  to  the  leading 
order  of  the  Chapman-Enskog  expansion,  and  one  considers  some  simple 
types  of  average  flows,  finding  an  average  flow  is  not  a  very  hard  task. 
For  instance,  in  the  practically  important  case  of  a  stationary  flow  in 

a  pipe  under  uniform  pressure  gradient,  it  is  easy  to  show  that,  u*  =  w' 

( i . e .  =  0  and,  hence.  dA  =  0  )  is  a  suitable  solution  around  which 

such  a  decomposition  can  be  carried  out. 

There  is  no  mathematically  sound  way  of  deciding  how  good  the  Chap¬ 
man  F.nskog  expansion  is  in  any  given  situation.  However,  intuitively 
one  expects  that  even  though  a  dilute  particle  phase,  if  left,  in  isola¬ 
tion.  would  be  quite  far  from  local  equilibrium,  an  interaction  with  the 
fluid  phase  will  bring  the  combined  system  quite  close  to  local  equili 
hrium  and  thus,  justify  the  use  of  a  truncated  Chapman-Enskog  expansion. 
To  some  extent,  this  would  justify  the  two-fluid  description  of  the  two 

5 

phase  flow;  set;  also  Drew  for  a  discussion.  It  is  to  be  emphasized, 
however,  that  our  approach  is  not  tied  to  the  use  of  such  an  expansion 

As  a  final  remark,  we  wish  to  draw  attention  to  the  amusing  fact, 
expressed  by  Kq .  (3.6),  namely  that  the  form  of  the  conservation  Jaws  is 
unchanged  despite  the  fact  that  the  velocity  dependence  of  the  forces 
invalidates  I.iouville's  theorem.  To  some  extent.,  this  situation  is 
similar  to  the  one  expressed  by  the  so  called  "Hilbert  paradox"  in 
transport  theory;  the  Boltzmann  distribution  function  contains  more 

information  than  the  totality  of  its  f ew  moments,  see  e.g.,  Uhlenheek.’ 


b 

D.  A.  Drew,  "Mathematical  Modelling  of  Two -Phase  Flow,"  in  Annua  1 
Review  of  Fluid  Mechanics,  Vo  1 .  15.  p  261.  Annual  Reviews,  Palo  Alto. 

1963 

f>6 .  K.  Uhlenheek,  "The  Validity  and  the  t,  imitations  of  the  Boltzmann 
Equation,"  Acta  Piiys ,  Austriaca.  Supp .  X,  pp  107  110,  1973. 
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This  Laboratory  undertakes  a  continuing  effort  to  improve  the  quality  of  the 
reports  it  publishes.  Your  comment s/answers  to  the  items/questions  below  will 
aid  us  in  our  efforts. 

1.  BRL  Report  Number _ _ _ Date  of  Report _ 

2.  Date  Report  Received _ _ .  . 

3.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or 

other  area  of  interest  for  which  the  report  will  be  used.) _ 


4.  How  specifically,  is  the  report  being  used?  (Information  source,  design 
data,  procedure,  source  of  ideas,  etc.) _ 


5.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far 
as  man-hours  or  dollars  saved,  operating  costs  avoided  or  efficiencies  achieved 
etc?  If  so,  please  elaborate. _ 


6.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future 
reports?  (Indicate  changes  to  organization,  technical  content,  format,  etc.) 
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